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B Meson Matrix Elements from Various Heavy Quark Effective Theories 

Terrence Draper * and Craig McNeile ^ 

^ Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA 

Various properties of heavy-light mesons are determined, including decay constants, the -Bu-parameter, and 
the Isgur-Wise function. The heavy (bottom) quark is simulated with the static, NRQCD and/or (fixed- velocity) 
lattice-HQET effective theories, using optimally-smeared sources as produced by the "Maximal Operator Smearing 
Technique" . 



1. INTRODUCTION 

Theoretical predictions of the B-meson decay 
constant (/s), the Bb parameter (from B^-B^ 
mixing), and the Isgur-Wise function ^ (the uni- 
versal form factor for semi-leptonic decays of 
heavy mesons), arc very important for narrowing 
the constraints on the CKM quark-mixing ma- 
trix, and must be determined by non-perturbative 
methods, of which the lattice regularization is the 
most promising. Current approaches include us- 
ing for the heavy quarks the same action (e.g. 
Wilson) as used for light quarks and hoping to 
control the lattice artifacts arising in the ex- 
trapolation to the regime Mga > 1, with a 
the lattice spacing. An alternative is the for- 
mulation of an effective theory for the heavy 
quarks directly on the lattice, such as the static 
theory pioneered by Eichten. The static ap- 
proximation and its non-zero velocity generaliza- 
tion, "lattice HQET" [1], used for the direct lat- 
tice calculation of the Isgur-Wise function, suffer 
from severely degraded signal-to-noise. To over- 
come this problem, we developed the variational 
technique "MOST" (Maximal Operator Smearing 
Technique) [2,3] which exploits all of the informa- 
tion available from relative smearing of the heavy 
and light quarks to construct optimal smearing 
functions. 
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Here we report that these same smearing func- 
tions ( "S-sources" ) , which allow for a determina- 
tion of /b in the static approximation with smah 
statistical errors and free from excited-state con- 
tamination [2,3], arc used successfully in a cal- 
culation of the Bfi-parameter (for static- Wilson 
mesons) , and for the calculation of the Isgur-Wise 
function. In addition, MOST can be used most 
effectively for the calculation of B meson prop- 
erties where the heavy quark is treated with the 
Non-Relativistic QCD effective theory. We report 
results for mass splittings and for decay constants 
using NRQCD through first order in l/Mg for 
both the action and the current [4,5]. All our re- 
sults are for 32 quenched /3 = 6.0 configurations 
on a 20^ x 30 lattice. 

2. STATIC Bb PARAMETER 

The four-Fermi operator contributing to B^- 

BO mixing is Ol = (67m(1 " 75)9)(&7m(1 " Ts)?)- 
Its matrix element normalized to its value in 
the vacuum-saturation approximation is the Bb- 
parameter 



SB(/i) 



(1) 



which would equal 1 if VSA were exact. The 
corresponding quantity in the lattice static effec- 
tive theory can be calculated on the lattice from 
a ratio of three- and two-point correlation func- 
tions for asymptotically large time separations. 
As shown in figure 1, the use of a MOST Z-source 
results in the early onset of a plateau with smah 
statistical errors. 
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Figure 1. Contribution to the Bb parameter from 
the ratio of three-point to two-point correlation 
functions which asymptotically equal the normal- 
ized matrix element of the Ol lattice operator. 

The continuum _BB-parameter is determined 
from the lattice static effective theory through 
a two-step matching procedure which results in 
mixing with new operators Os, On, and Or [6]. 
Lattice values for the raw B-parameters for each 
of these operators are calculated analogously to 
that of Ol, and are tabulated in table 1. The 
matching coefficients have been calculated for 
Bsfl by Flynn, Hernandez and Hill [6]. We 
use their expressions with the following revisions: 
we subtract perturbative corrections for to 
get matching coefficients for B itself, we use 
the Lepage-MacKenzie prescription [7] to get a 
"boosted" value for the gauge coupling evaluated 
at the appropriate momentum scale of the pro- 
cess, and we use renormalization- group-improved 
perturbation theory [8] . Our preliminary result is 
BBimt,) = imtl (Sij(2.0GeV) = 1.12l^.) For 
details of the calculation, see [9]. 

3. ISGUR-WISE FUNCTION 

We calculate the Isgur-Wise function using 
the lattice HQET formalism of Mandula and 
Ogilvie [1] which seeks to formulate the heavy- 
quark effective theory directly on the lattice. 
Their initial attempts [1] suffered from a lack of a 
convincing signal, as the signal-to-noise degraded 
before ground-state saturation. Since at zero ve- 
locity the lattice HQET reduces to the static the- 
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Table 1 

Raw lattice B parameters for the various opera- 
tors appearing in the continuum-lattice matching, 
and the B parameter for the linear combination 
of operators corresponding to the continuum op- 
erator O. 



ory, for which we found that ''MOST' was very 
effective in alleviating the problem of ground- 
state extraction, it is natural to apply MOST to 
HQET. In fact, for low velocities, which is what 
is wanted to extract the slope of the Isgur-Wise 
function near the normalization point, it suffices 
to use as sources the smeared sources pro- 
duced by MOST for zero velocity. 

This is demonstrated in figure 2 for which the 
"effective mass" has plateaued at very early time 
separations, so that a reliable ground state signal 
is seen well before the signal-to-noise ratio has 
degraded. 

Following Mandula and Ogilvie [1], the (UN- 
renormalized) Isgur-Wise form factor is extracted 
from a ratio of three-point functions. From fig- 
ures 3, 4 and table 2, it appears that MOST's Z- 
sources are much better at removing excited state 
contaminations than previous choices and that, 
for the first time, a result with about 5% statisti- 
cal errors and 5% systematic errors (from excited- 
state contamination only) is available. We do not 
quote a physical quantity, however, since our re- 
sult is not yet renormalized. The quality of the 
present numerical calculation justifies proceeding 
to a calculation of the renormalization which is 
forthcoming. 
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Figure 2. Effective mass plot for the (velocity 
V = 0.25) HQET-light meson two-point correla- 
tion function (local sink, smeared (rest frame) Z- 
source). 
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Table 2 

Negative of the slope at the normaHzation point, 
— ^ (1), from the ratio of three-point functions 
which equals the (UN-renormalized) Isgur-Wise 
function (^{v-v ) at asymptotically-large times At. 



Figure 3. (UN-renormalized) Isgur-Wise function 
for t2 — ti = ti — to = At = 2 (local sink, smeared 
(rest frame) .Z-source). 
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Figure 4. Same as for figure 3 except for At = 3. 



4. NRQCD 

With the success of NRQCD in describing 
heavy quarkonia, it is natural to explore its use 
for heavy-light systems. Here we report results on 
decay constants and mass splittings for i?-mesons 
with NRQCD heavy quarks (tadpole improved) 
and Wilson light quarks in the quenched approx- 
imation. We upgrade earlier calculations [3], by 
including terms through 1/Mq in the axial cur- 
rent as well as in the Lagrangian, as was done 
in the calculations [4] for NRQCD- Wilson with 
dynamical staggered quarks, and for quenched 
NRQCD-clover [5]. 



MOST is not only feasible but very effective for 
NRQCD-light mesons, as is demonstrated in fig- 
ure 5, where again we see an early plateau in the 
effective mass plot. Table 3 lists results for the 
iriB^-iTiB splitting which agrees with our earlier 
calculation in the static approxiniation (as pre- 
dicted by HQET); the splitting is smaller than 
experiment (96(6) MeV). The mB.-ms splitting 
is smaller than experiment (46(1) MeV) and than 
with a clover Hght quark (37(6) MeV) [5]. 

Figure 6 shows a ratio of two-point corre- 
lation functions which should approach = 
(/b-\/™b/2)'^** for large times. We see a plateau 
immediately, which again indicates the efficacy of 
MOST in isolating the ground state. 
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Table 4 _ 

I I'l (UN-renormalized) = {fB^/rriBj^Y^^^ for Qs and Qd pscudoscalar mesons with Wilson light 
quarks and NRQCD heavy quarks ( with tadpole improvement). The presence (absence) of the super- 
script [°1 indicates that the 0(1/M) correction is not (is) included in the matrix element of the current. 

APE-Iitude for 0(1/M) corrected current 



"tl=0" 
"tl=l" 
"tl=2" 
"tl=3" 
"tl=4" 

i I I i j i i i i i 



20"3 ' 30, beta=6.0, 32 ctg. 

NRQCD: (coUab.) K.E. + sigma.B, synim. update 

11=2, Ma=1.7, u_0=0.88 

Kappa = 0.154 

Z[a#=3;t_0=2,t=3] 





mBs —iriB 


mpt. -TUB 


(2.1 GeV ) 


(2.1 GeV ) 


(static) 00 


78I4 MeV 





2.0 


76l^ MeV 


30l^ MeV 


1.7 


771^ MeV 


35l^ MeV 



Table 3 

Bg B and B* B mass splittings. 
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Figure 5. Effective mass plot for the pseudoscalar. 

But this ratio uses unnecessarily large times, 
and our results, presented in table 4, are from 
simultaneous fits to local-smeared and smeared- 
smcarcd two-point correlation functions. Our ra- 
tio fsj fs = 1.151? compares with the following 
quenched /3 = 6.0 results: 1.18(1) [10] (extrapo- 
lated Wilson heavy quark) and 1.12 [5] (NRQCD 
heavy and clover hght). We got 1.221J [3] for 
static-heavy-Wilson-light (without tadpole im- 
provement) . 



Time (t2) 

Figure 6. "APE" ratio of LS and SS correlation 
functions for the axial current. 

REFERENCES 

1. J. Mandula and M. Ogilvie, Nucl. Phys. B 
(Proc. Suppl. ) 34, 480 (1994). 

2. T. Draper and C. McNeile, JVucJ. Pbys. B 
(Proc. Suppl. ) 34, 453 (1994). 

3. T. Draper, C. McNeile and C. Nenkov, JVucJ. 
Phys. B (Proc. Suppl. ) 42, 325 (1995). 

4. S. CoUins et al, Nucl. Pbys. B (Proc. Suppl. 
) 42, 395 (1995). 

5. A. Ali Khan et al, Nucl. Phys. B (Proc. 
Suppl. ) 42, 382 (1995). 

6. J. Flynn, O. Hernandez and B. Hill, Phys. 
Rev. D 43, 3709 (1991). 

7. P. Lepage and P. Mackenzie, Phys. Rev. D 
48, 2250 (1993). 

8. V. Gimenez, JVucJ. Phys. B401, 116 (1993). 

9. J. Christensen, T. Draper and C. McNeile, in 
preparation. 

10. (MILC cohaboration) C. Bernard et al, Nucl. 
Phys. B (Proc. Suppl. ) 42, 388 (1995). 



